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DESCRIPTION 

METHOD AND APPARATUS FOR PREDICTING BENDING LIFE SPANS OF 
ELECTRIC WIRES AND/OR WIRE PROTECTING MEMBERS INDUCED BY 
VIBRATIONS, AND RECORDING MEDIUM STORING PROGRAM 

5 

Technical Field 

The present invention relates to a method and apparatus 
for predicting bending life spans of a plurality of electric 
wires clamped at least two points to a predetermined part and/or 
10 their wire protecting members, the bending induced by vibrations, 
and a program for predicting the same. 

Background Art 

Usually, in a vehicle or the like, a plurality of electric 

15 devices are electrically connected to one another by a plurality 
of electric wires (referred to frequently and simply as wires) , 
The wires are bundled, by an insulation locking band, tape or 
the like, into a single bundle as so-called wire harness which 
is extendedly arranged at predetermined locations of the 

20 automobile or the like. In some cases, the wires are protected 
by a wire protecting member such as a grommet are arranged at 
a predetermined part. In other cases, the wires are arranged 
at a predetermined part of the vehicle or the like, while not 
bundled into one line and/or not protected by the wire protecting 

25 members. 
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The vehicle is always placed in such an environment that 
it is always vibrated due to engine drive and the like. In 
particular, the electric wires distributed inside the engine 
room directly receiving the engine vibrations, are repeatedly 
5 bent and deformed, and finally disconnect due to the vibrations . 
Experimental results of such are reported. In this respect, 
it is of particular importance to accurately predict the bending 
life spans of the electric wires and the wire protecting members . 

Generally, the prediction of the bending life spans of 

10 the electric wires and the wire protecting members is carried 
out by repeating design, trial manufacturing, and endurance 
test. For example, in such an endurance test, wires distributed 
along an envisioned wiring path are placed on a vibration base 
plate . The vibration base plate is vibrated at a predetermined 

15 frequency and a predetermined amplitude by the vibrator, and 
a bending life span of the wire when it is vibrated at a specific 
number of vibrations is predicted. 

The literature cited in the present specification are: 
Non-patent document 1 

20 "Matrix Infinite Element Method" by B. Nasu, published 

by Brain Library Publishing Corp, August 10, 1978, pp 7 to 15 
Non-patent document 2 

"Mode Analysis and Dynamic Design" by Kimihiko Yasuda, 
published by Corona Corp, November 10, 1993, pp 54 to 56 

25 On the other hand, a development period of the vehicle 



tends to be shorter, and further there is a demand of enhancing 
a bending life prediction accuracy- In this circumstance, the 
related bending life span predicting method in which the 
endurance test is repeated cannot satisfactorily reduce the 
5 development period, and cannot sufficiently meet the demand 
of enhancing the prediction accuracy. 

Disclosure of the Invention 

Accordingly, to cope with the above-mentioned problems, 
10 an object of the present invention is to provide a method and 
apparatus for predicting bending life spans of electric wires, 
which can meet the demands of a highly accurate prediction of 
the bending life span and of reducing the development period 
in an environment where vibrations occur, and a recording medium 
15 storing a program for predicting the same in such ways. 

In order to solve the aforesaid object, the invention 
is characterized by having the following arrangement. 

(1) A bending life predicting method of predicting bending 
life span of a plurality of wires induced by vibration, at least 
20 two points of each of the plurality of wires being constrained, 
the method comprising the steps of: 

a pre-storing step of pre-storing each predicting 
function representing relationships among atmosphere 
temperatures, stresses and bending life spans for the plurality 
25 of wires; 




a setting step of setting the plurality of wires, the 
atmosphere temperatures, pre-vibration shapes of the plurality 
of wires, and constraint conditions of the plurality of wires; 

a finite element model forming step of forming finite 
5 element models of the plurality of wires by using a finite element 
method; 

a vibration analyzing step of calculating natural 
frequencies for the pre-vibration shapes and calculating 
stresses in individual finite elements of the finite element 
10 models which correspond to the natural frequencies, 
respectively; 

a maximum stress retrieving step of retrieving a maximum 
stress from the stresses calculated in the vibration analyzing 
step, for each of the plurality of wires; 
15 a predicting function readout step of reading predicting 

functions corresponding to the atmosphere temperatures set in 
the setting step, respectively; 

a bending life predicting step of acquiring a bending 
life span corresponding to the maximum stress of each of the 
20 plurality of wires while referring to the predicting functions 
read out in the predicting function readout step, and obtaining 
a shortest bending life span from the bending life spans; and 

an output step of outputting the shortest bending life 
span obtained in the bending life predicting step. 

25 
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(2) The method according to (1), wherein in the vibration 
analyzing step, the plurality of wires are regarded as a wiring 
structure in which the plurality of wires are bundled, and 
natural frequencies of the plurality of wires are computed, 

5 respectively. 

(3) The method according to (1), wherein the plurality of 
wires are bundled into a single bundle, and the bundling of 
the plurality of wires is set as one of the constraint conditions 

10 in the setting step. 

(4) The method according to (1), further comprising aposition 
specifying step of specifying a position on the wire 
corresponding to the shortest bending life span, the output 

15 step outputting the position specified by the position 
specifying step. 

(5) The method according to (1), wherein 

in the vibration analyzing step, displacements of finite 
20 elements of the finite element models which correspond to the 
natural frequencies are calculated, 

the method includes an interference part predicting step 
of predicting an interference part on the plurality of wires 
which is induced by vibrations based on the calculated 
25 displacements, and 




the output step outputs the predicted interference part. 

(6) The method according to (1), wherein a curve representing 
a lower confidence interval to a population regression function 
5 statistically calculated using the stresses and data on bending 
endurance life spans that are obtained under a plurality of 
typical atmosphere temperatures for the plurality of wires, 
is used. for the predicting function. 

10 (7) A bending life predicting method of predicting bending 
life spans of a plurality of electric wires and wire protecting 
members for protecting the plurality of wires from being bent 
induced by vibrations, at least two points of each of the 
plurality of wires being constrained, the method comprising: 

15 a pre-storing step of pre-storing relationships among 

atmosphere temperatures, stresses and bending life spans for 
plurality of wires and the wire protecting member; 

a setting step of setting the plurality of wires, the 
wire protecting member, the atmosphere temperatures, 

20 pre-vibration shapes of the plurality of wires and the wire 
protecting member, and constraint conditions of the plurality 
of wires and the wire protecting member; 

a finite element model forming step of forming finite 
element models of the plurality of wires and the wire protecting 

25 member by using a finite element method; 




a vibration analyzing step of calculating natural 
frequencies for the pre-vibration shapes of the plurality of 
wire and the wire protecting member, and calculating stresses 
in individual finite elements of the finite element models which 
5 correspond to the natural frequencies, respectively; 

a maximum stress retrieving step of retrieving maximum 
stresses from the stresses calculated in the vibration analyzing 
step, for each of the plurality of wires and the wire protecting 
member; 

10 a predicting function readout step of reading predicting 

functions corresponding to atmosphere temperatures set in the 
setting step; 

a bending life predicting step of acquiring bending life 
spans corresponding to the maximum stresses of the plurality 
15 of wires and the wire protecting member, while referring to 
the predicting functions read out in the predicting function 
readout step, respectively, and obtaining a shortest bending 
life span from said bending life spans; and 

an output step of outputting the shortest bending life 
20 span obtained in the bending life predicting step. 

(8) The method according to (7), wherein in the vibration 

analyzing step, the plurality of wires are regarded as a wiring 

structure in which the' plurality of wires are bundled, and 

25 natural frequencies of the plurality of wires are computed, 
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respectively. 

(9) The method according to (7), wherein the plurality of 
wires are bundled into a single bundle, and the bundling of 

5 the plurality of wires is set as one of the constraint conditions 
in the setting step. 

(10) The method according to (7), further comprising a position 
specifying step of specifying a position on the wire or the 

10 wire protecting member corresponding to the shortest bending 
life span, the output step outputting the position specified 
by the position specifying step. 

(11) The method according to (7), wherein 

15 in the vibration analyzing step, displacements of finite 

elements of the finite element models which correspond to the 
natural frequencies are calculated, 

the method includes an interference part predicting step 
of predicting an interference part on the plurality of wires 
20 which or the wire protecting member is induced by vibrations 
based on the calculated displacements, and 

the output step outputs the predicted interference part. 



25 



(12) The method according to (7), wherein a curve representing 

a lower confidence interval to a population regression function 
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statistically calculated using the stresses and data on bending 
endurance life spans that are obtained under a plurality of 
typical atmosphere temperatures for the plurality of wires and 
the wire protecting member, is used for the predicting function. 

5 

(13) A bending life predicting device for predicting bending 
life spans of a plurality of wires induced by vibrations, at 
least two points of each of the plurality of wires being 
constrained, the device comprising: 
10 a pre-storing unit for pre-storing each predicting 

function representing relationships among atmosphere 
temperatures, stresses and bending life spans for the plurality 
of wires; 

a setting unit for setting the plurality for wires, the 
15 atmosphere temperatures, pre-vibration shapes of the plurality 
of wires, and constraint conditions of the plurality of wires; 

a finite element model forming unit for forming finite 
element models of theplurality of wiresbyusing a finite element 
method; 

20 a vibration analyzing unit for calculating natural 

frequencies for the pre-vibration shapes and calculating 
stresses in individual finite elements of the finite element 
models which correspond to the natural frequencies, 
respectively; 

25 a maximum stress retrieving unit for retrieving a maximum 




stress from the stresses calculated in the vibration analyzing 
unit, for each of the plurality of wires; 

a predicting function readout unit for reading predicting 
functions corresponding to the atmosphere temperatures set in 
5 the setting unit, respectively; 

a bending life predicting unit for acquiring a bending 
life span corresponding to the maximum stress of each of the 
plurality of wires while referring to the predicting functions 
re.ad out in the predicting function readout unit, and obtaining 
10 a shortest bending life span from the bending life spans; and 

an output unit for outputting the shortest bending life 
span obtained in the bending life predicting unit, 

(14) The device according to (13), wherein 

15 the vibration analyzing unit calculates displacements 

of finite elements of the finite element models which correspond 
to the natural frequencies, 

the device includes an interference part predicting unit 
for predicting an interference part on the plurality of wires 
20 which or the wire protecting member is induced by vibrations 
based on the calculated displacements, and 

the output unit outputs the predicted interference part. 

(15) A computer readable recording medium storing a program 

25 f or predictingbending life spans of a plurality of wires induced 
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by vibrations, at least two points of each of the plurality 
of wires being constrained, the program causing a computer to 
function as: 

a pre-storing unit for pre-storing each predicting 
5 function representing relationships among atmosphere 
temperatures, stresses andbending lif e spans for the plurality 
of wires; 

a setting unit for setting the plurality for wires, the 
atmosphere temperatures, pre-vibration shapes of the plurality 
10 of wires, and constraint conditions of the plurality of wires; 

a finite element model forming unit for forming finite 
element modelsof the plurality of wires by using a finite element 
method; 

a vibration analyzing unit for calculating natural 
15 frequencies for the pre-vibration shapes and calculating 
stresses in individual finite elements of the finite element 
models which correspond to the natural frequencies, 
respectively; 

a maximum stress retrieving unit for retrieving a maximum 
20 stress from the stresses calculated in the vibration analyzing 
unit, for each of the plurality of wires; 

a predicting function readout unit for readingpredicting 
functions corresponding to the atmosphere temperatures set in 
the setting unit, respectively; 
25 a bending life predicting unit for acquiring a bending 



life span corresponding to the maximum stress of each of the 
plurality of wires while referring to the predicting functions 
read out in the predicting function readout unit, and obtaining 
a shortest bending life span from the bending life spans; and 
5 an output unit for outputting the shortest bending life 

span obtained in the bending life predicting unit. 

Brief Description of the Drawings 

Fig. 1 is a block diagram showing a basic configuration 
10 of the present invention. 

Fig. 2A is a diagram showing an electric wire handled 
in the invention when it is discretized. Fig. 2B is a diagram 
showing beam elements and nodes of the discrete electric wire 
shown in Fig. 2A. 
15 Fig. 3 is a block diagram showing a hardware configuration 

constructed according to the invention. 

Figs. 4A and 4B are diagrams showing bending life span 
data stored in a storing device in Fig. 3. 

Fig . 5 is a flow chart showing a main processing procedure 
20 according to an embodiment of the invention. 

Fig. 6 is a flow chart showing a vibration analysis, 
processing procedure shown in Fig. 5. 

Figs . 7A, 7B and 7C are diagrams showing predictive shapes 
of the electricwire inprimary, secondary and tertiary vibration 
25 modes, respectively. 
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Best Mode for Carrying out the Invention 

The preferred embodiment according to the present 
invention will be described with reference to the accompanying 
5 drawings. Pre-conditions, theories and basic formulae which 
are used in the embodiment will first be described with reference 
to Figs. 2A and 2B. Fig. 2A is a diagram showing an electric 
wire in a state it is discretized, handled in the invention. 
Fig. 2B is a diagram showing beam elements and nodes of the 
10 discrete electric wire shown in Fig. 2A. 

In the invention, a finite element method is utilized 
on the basis of the following assumptions: 

(1) The electric wire consists of an elastic member. 

(2) The electric wire is an aggregation of beam elements 
15 serially connected to one another. 

(3) A linearity is retained in each beam element. 

(4) A cross section of the electric wire is circular. 

Next, this electricwire 1 is discretized. Specifically, 

20 as shown in Fig. 2A, the electric wire 1 is divided (discretized) 

into plural beam elements CI, C2, C3, ... . In other words, the 

electric wire 1 is regarded as an aggregation of a finite number 

of beam elements serially connected to one another. 

Accordingly, the electric wire 1 may be considered as 

25 an aggregation of beam elements CI, C2, C3, ... which are serially 
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connected at nodes Nl, N2, N3, ... to one another, as shown in 
Fig. 2B. Properties necessary for the beam element are, for 
example, length "1" (see Fig. 2A) , a cross section area A (see 
Fig. 2A) , a moment of inertia of area, a polar moment of inertia 
5 of area, density, amodulus of longitudinal elasticity, amodulus 
of transverse elasticity, and the like. Those properties may 
be obtained in advance by measurement or calculation. It is 
preferable that the properties are formed as a database, and 
available ad libitum. In the specification, the length and 

10 the cross sectional area are set as geometric properties, and 
the moment of inertia of area, the polar moment of inertia of 
area, the modulus of longitudinal elasticity, and the modulus 
of transverse elasticity are set as material properties. 

As well known. Hook's Law holds for a displacement of 

15 a structure not being vibrating by a static force within a range 
of elasticity of the structure. Likewise, Hook's Law holds 
also for the beam elements CI, C2, C3, as is known. The 
conditions of the displacement continuity and the equilibrium 
of forces are satisfied among the adjacent beam elements CI, 

20 C2, C3, .... By utilizing those facts, those beam elements CI, 
C2, C3, ... are connected and the resultant discrete body (electric 
wire) is expressed by the following expression (1) . 



[K] {X} = {F} ... (1) 

25 In the above expression, [K] is a global stiffness matrix which 
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are assembled from stiffness matrices of individual beam 
elements. To state it simply, the contents of the global 
stiffness matrix [K] is like an aggregation of the properties. 
Those values are previously set. {x} is a nodal displacement 
5 vector (or a displacement vector, simply) , and an array of 
displacement components of all the nodes. The component of 
the nodal displacement vector {x} includes a preset vector, 
such as a point of constraint, and an unknown quantity. {F} 
is a net external force vector, and an array of external force 

10 components at all the nodes. 

In this way, a finite element model of the discrete wire 
having divided into an optional number of beam elements is formed. 
A technique, which is similar to a technique disclosed in 
Japanese patent applicationNo, 2002-279503, filed on September 

15 25, 2002 by the Applicant of the present patent application, 
can be applied to the formation of the finite element model. 
The general matrix finite element method is described in, for 
example, the non-patent documents described above. 

In the invention, to predict a bending life span of the 

20 electric wire induced by vibrations, it is necessary to obtain 
the natural frequencies, displacements, stress and the like. 
A theory and calculation formulae, which are used for obtaining 
the natural frequencies, displacements, stress and the like, 
will be brief lydescribedbelow, while developing the expression 

25 (1). 




To dynamically handle the wire for predicting its bending 
life span, an inertia force [M] {x"} is added to the load term 
in the expression of the static equilibrium expressed by the 
expression (1), and damping [C] {x' } is taken into account, 
5 and then we have 

[M] {X"} + [C] {X'} + [K] {X} = {F (t)} ... (2) 

where 

[M] : inertia matrix 
10 [C] : damping coefficient matrix 

{x' } : first order differential of the displacement vector {x} 
{x"} : second order differential of the displacement vector 

{X} 

15 In a free vibration state, if external force = 0 ({F(t) } 

= 0 in the expression (2) ) , and no damping is present [C] {x' } 
= 0 in the expression (2) ) , each point of the wire oscillates 
in a simple harmonic motion at a fixed frequency. When the 
following expression (3) is substituted into the expression 

20 (2), 

{x} = {Xo} sin Qt ... (3) 

then we have 

[M] ({Xo} sin Ot)" + [K] ({xo} sin Qt) = 0 
(-co^ [M] + [K] ) {x} sing cot = 0 
25 ( [K] - CD^ [M]) {X} = 0 ... (4) 




From the expression (4), the natural frequency co and the 
displacement vector {x} as an eigen-vector can be obtained. 
Even if the displacement vector {x} is multiplied by a fixed 
5 number, the expression (4) holds. Hence, the displacement 
vector {x} is a relative displacement. Accordingly, a stress 
obtained from the displacement vector {x} is also a relative 
value. An absolute value obtained from the expression (4) is 

only the natural frequency co . To calculate actual 

10 displacements and stresses, it is necessary to input a load 
condition of an actual vibration. 

To dynamically handle the wire for its bending life span 
prediction, if the time-dependant load {F (t) } is input to the 
expression (2) , displacements and stresses each for each unit 

15 time can be calculated, basically. In case where the equation 
is solved by inputting the time-dependant load {F (t) } to the 
expression (2), an analysis time to reproduce the vibration 
is long in a phenomenon attendant with periodic vibrations, 
and hence, a time to output the analysis result is extremely 

20 long. To avoid this, a common practice is to convert the 
expression from a time domain to a frequency domain by Fourier 
transformation. In the calculation of the expression in the 
frequency domain, viz., in the frequency response analysis, 
the frequency-dependant load is input to the expression to 

25 calculate displacement and stress per each frequency. For the 




frequency-dependant load, a force applied to points of 
constraint on a wire harness or a wire protecting member as 
an object under prediction of the bending life is set at an 
optional value, the force attending envisioned vibrations of 
5 the engine, vehicle body, and the like. Thereafter, the 
expression is returned from the frequency domain to the time 
domain by inverse Fourier transformation, whereby an amount 
of calculation is reduced. This technique is well known, and 
hence, no further description of it will be given here. 

10 By utilizing the theories and calculation formula thus 

far described, one can calculate natural frequency, 
displacements and stresses, which will be described later on. 
The general eigenvalue analysis as described above is described 
also in the non-patent document 2, for example. 

15 In the present invention, a vibration analysis of a 

plurality of electric wires including, for example, grommets 
as wire protecting members is also carried out. To apply the 
finite element method to the wire protecting member, the wire 
protecting member is divided into a plurality of triangular 

20 finite elements, and stresses at those finite elements are 
obtained. A vibration analysis of the wire protecting member, 
which follows the formation of the finite element model 
consisting of such finite elements, is carried out in a manner 
as in the case of the electric wire. 

25 Next, a hardware configuration for realizing a process 




procedure constructed according to the present invention will 
be described. Fig. 3 is a block diagram showing a hardware 
configuration constructed according to the invention. Figs. 
4A and 4B are diagrams showing bending life span data stored 
5 in a storing device in Fig. 3. 

As shown in Fig. 3, in the present invention, a personal 
computer, for example, is used which has a basic configuration 
made up of a micro-computer 51, an input device 52, a display 
device 53, a printing device 54, a storing device 55, a 

10 communication interface 56, and a read/write device 57. The 
micro-computer 51 includes a CPU (central processing unit) 51a, 
a ROM 51b for storing a boot program and others, and a RAM 51c 
for temporarily storing various processing results . The input 
device 52 is a keyboard, a computer mouse and the like; the 

15 display device 53 is an LCD, a CRT and the like for displaying 
processing results; and the printing device 54 is a printer 
for printing processing results. 

The storing device 55 is, for example, a hard disc drive, 
and the communication interface 56 is, for example, a modem 

20 board for communicating with external devices by way of an 

internet, a L7\N line and the like. The read/write device 57 

reads a bending life span predicting program 59a constructed 

according to the invention, which is stored in a recording medium 

59, and writes a result file 55c into the recording medium 59. 

25 Those constituent components are interconnected through an 
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internal bus 58. 

At least a set data file 55a, a bending life data file 
55b, and a result file 55c are stored in the storing device 
55. The set data file 55a contains data preset for predicting 
5 the bending life span. The data represents types of electric 
wires and wire protecting meinbers, atmosphere temperatures, 
pre-vibration shapes, constraint conditions and others, which 
are set in a step SI to be described later. The set data file 
55a contains position information of obstacles 2, such as 

10 pre-estimated stays and electric parts (see Fig. 7) . 

The bending life data file 55b is an aggregation of 
predicting functions yl, y2 and y3 as statistically calculated 
using data of stresses and the number of bendings for endurance 
on electric wires 55al and 55a2, and wire protecting member 

15 (e.g., grommet) 55a3, which are gathered at different atmosphere 
temperatures of, for example, - 40°C, 0°C and 25 °C. 

A curve y22 represents an upper confidence interval and 
the curve y23 represents a lower confidence interval to a 
population regression function y21 obtained by the known 

20 regression analysis, as shown in Fig. 4B. For the predicting 
function, the curve y23 of those curves y22 and y23 is used. 
The confidence interval is, for example, 95%. Such predicting 
functions are obtained under condition that the electric wire 
and the grommet are placed at different atmosphere temperatures . 

25 Accordingly, the bending life span is predicted under more severe 
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statistical conditions. As a matter of course, a given 
statistical reliability on the predicting function is taken 
into account, and a computing process of the predicting function 
iseasy. Asaresult, thebending life span is severelypredicted 
5 without addition of complicated processing procedures. This 
contributes to further enhancement of quality and wiring path 
plan. Incidentally, the number of bendings for endurance may 
be obtained using the population aggression function. 

The result file 55c contains all the stresses of the 
10 discrete wire and the discrete grommet in the finite elements, 
which are recorded by each natural frequency. The result file 
55c is stored in a text format, and may be output ad libitum. 
The storing device 55 corresponds topredicting function storing 
means set forth in claims. 

With such a configuration, the micro- computer 51 installs 
a bending life span predicting program 59a that is read by the 
read/write device 57 into the storing device 55. After power 
on, the micro-computer 51 is activated according to a boot 
program that is stored in the ROM 51b, and starts the bending 
life span predicting program 59a. According to the bending 
life span predictingprogram 59a, themicro-computer 51 predicts 
a bending life span of the electric wire and/or the wire 
protecting member induced by vibrations, causes the display 
device 53 and the printing device 54 to display and print the 
prediction result, and stores the result in the storing device 
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55. The bending life span predicting program 59a may be 
installed to another personal computer having the same 
configuration, and after installing, operates the computer as 
a bending life predicting device. The bending life span 
5 predicting program 59a may be provided through a communication 
line, such as an Internet, a LAN or the like, instead of the 
recording medium 59. 

Processing procedures according to the present 
embodiment of the invention will be described using an 
10 explanatory diagram of Figs. 7A through 7C and flow charts of 
Figs . 5 and 6 . Fig . 5 is a flow chart showing a main processing 
procedure according to an embodiment of the invention. Fig. 
6 is a flow chart showing a vibration analysis processing 
procedure Shown in Fig. 5. Figs. 7A, 7B and 7C are diagrams 
15 Showing predictive shapes of the electric wire in primary, 
secondary and tertiary vibration modes, respectively. 

As shown in Fig. 5, data necessary for the prediction 
is set in a step SI . Specifically, items to be set are at least 
a plurality of electric wires (exactly types of electric wires) 
whose bending life spans are to be predicted, wire protecting 
members, (sometimes not required), atmosphere temperatures, 
pre-vibration shapes for each wires, and constraint conditions 
for each electric wires . To set those necessary items, an input 
screen (not shown) is presented on the display device 53, and 
25 necessary items are entered to the input screen by use of the 
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input device 52. 

The types of the electric wires are th« ^ 

wires are the types of the 

electric wires whose bending life spans ar. i- k 

y -Li-e spans are to be predicted 

For the electric wires, geometrical properti.. . 

properties and material 

properties ere Mn.ed to the type. o. t.e electric „ire= 
Speci.icaaa,, the ,eo™etric properties inciuCe indorsation 
about lengths an. cross sections .or each wire. The material 
properties are a „o»ent o. inertia o. area, a polar .o„ent o. 
rnertia of area, Censity a so.uius o. ion,it„.inaX elasticity 
10 an. a so.ui.s o. transverse elasticity. Those ^teriais „a 

be obtained in advance by tests for- 

oy tests, for example. Those properties 

concerns the finite elements in the stifi^n 

cne stiffness matrix [K] of 

the expression a,. The atmosphere temperatures are 
temperatures aroun. the eiectric wires an./or wire protecting 
.embers when the hendin, ii.e span prediction is carried out, 
and typically -40°c, 0°c and 250c. 

The pre-vihration shape is a shape of the electric wire 
or the li.e in a static state as indicated by reference numeral 
-inKi..,.. -theillustration,itissimpli.iedanddepicted 
irnearly. Pre-vihration shape la may .e manually entered 

by use Of the input device 52, or a pre-vi.ration shape that 
was previously forced before the bendin, life span prediction 
be used for the sa.e. Preferably, such a pre-vi.ration 

Shape as to satisfy given con^i-7-=.-i r,^- 

^ constraint conditions and the minimum 
bending radius of the elpr-i-^i^ . • 

the electric wire is calculated in advance, 
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and it is used for the pre-vibration shape la. 

For the constraint conditions, information of perfect 
constraint, rotation constraint and perfect free are set, while 
corresponding to coordinate information on nodes nl to n5 in 
5 fig 7A. The nodes nl to n5 are connection points of the finite 
elements into which the electric wire is divided. Some of the 

nodes correspondinpositiontopartsatwhichconstraintmembers, 
such as connectors and clips, respectively, m this instance, 
the nodes nl and n5 are put under the perfect or rotation 
10 constraint condition by the constraint members, while the 
remaining nodes n2 ton4 are put under the perfect free condition . 
The values set here concern the values of the finite elements 
in the displacement vector {x} in the expression (1) . 

in some cases, the plural electric wires which form a 
15 wire harness are bundled, by an insulation locking band, tape 
or the like, into a single bundle, and in other cases, those 
electric wires are not bundled, m case where the electric 
wires are bundled together, the wire bundling is set as the 
constraint condition, m setting the constraint condition of 
20 the wire bundling, it is preferable to handle the electric wires 
that come in contact with the insulation locking band or tape 
in distinction from those located on the inner side of the former 
electric wires. By so doing, the bending life spans of the 
Wires can be more exactly predicted since the electric wires 
25 bundled and the electric wires not bundled are separately 
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in some cases, those electric wires are passed through 
a Wire protecting „en*er, such as a gro,»et, and in other cases, 
the electric wires are arranged or distributed without using 
the Wire protecting meinher. The plural electric wires passed 
through the wire protecting raeinber are used such that as 
described above, in some instances, those wires are bundled 
together into a single bundle, and in other instances, those 
"ires are not bundled. In case where the wire protecting member 
is used, a constraint condition suitable for the wire protecting 
mender is to be set . The step SI thus far described corresponds 
to "setting step- and "setting unit" in claims. 

Next, in a step S2, finite element models of the electric 
Wires and/or wire protecting members are formed. The finite 
element model of the wire protectingmember is formed by dividing 
the wire protecting member into plural rectangular finite 
elements. The finite element model of the electric wire is 
formed such that the electric wire is divided into plural beam 
elements. The step S2 corresponds to a finite element model 
forming step and a finite element model forming means in claims . 

Next, m a step S3, a vibration analyzing process is 
carried out for each finite element. Specifically, m a step 
S301, on the assumption that a wire harness as a bundle of a 
Plurality of electric wires as objects vibrates, the expression 
W applies to the wire harness to calculate a natural frequency 
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of the wire harness, in this ea<;o 

n cms case, the geometric properties 

and the material properties as the individual elements of the 
stiffnessmatrix (K, calculated on the assumption that those 
properties are the results of composing the geometricproperties 
and the material properties of each electric wire and/or wire 

protecting member. For the nai-ii^=T ^ 

Zhe natural frequency, its values 

corresponding to, for example, thefir-.,i- 

P-Le, tne first, second and tertiary 

Vibration modes are calculated as shown in Figs . 7A to 7c The 
number of vibration modes to be calculated and the natural 
frequencies are not limited to the illustrated ones. Por 
example, the first, second and tertiary vibration modes are 
set as default values, and may be changed by use of the input 
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Following the natural frequency calculation, steps S302 
and S303 are executed in which the natural frequencies are 
applied to an expression obtained by Fourier transforming the 
expression ,2) as described above to thereby calculate 

displacements and stresses of all the f inn 1-0 oi 

dxj. tne finite elements for each 

electric wire and/or wire protecting member, m this case 
information as to whether or not the electric wires are bundled 
.nto one bundle and information as to whether or not specific 
electric wire or wires come in contact with the insulation 

locking band or the tape, which are set as th» . ■ 

i^eu as the constraint 

conditions in the <^tc^r^ qi 

the step SI, are incorporated regarding as the 

informationconcernallof the inertia! ..=4- • r 

^ne inertial matrix [MJ , the damping 




coefficient matrix [c], and the external force {F (t) ) in the 
expression (2) . 

In a step S304, the displacements and the stresses 
calculated every each natural frequency are linked with 
5 (finite elements) information of the position of each electric 
wire and/or wire protecting member as an object under prediction 
of thebendinglife, and are written into the result file . Those 
steps S3, S301 to S303 correspond to a vibration analyzing step 
and vibration analyzing unit in claims. 
10 subsequently, in a step S4, a maximum stress is retrieved 

from those stresses recorded in the result file, for each 
electric wire and/or wire protecting .uember under prediction 
of the bending life. In a step S5, a predicting function to 

theelectricwireand/orwireprotectlngmemberunderprediction 
Of the bending life is read out. Specifically, a predicting 
function of the atmosphere temperature as set in the step SI 
to the electric wire and/or wire protecting member under 
prediction of the bending life is read out from the bending 
life data file 55b stored in the storing device 55. The step 
S4 corresponds to a maximum stress retrieving step and a maximum 
stress retrieving unit in claims, and the step S5 corresponds 
to a predicting function readout step and predicting function 
readout unit in claims. 

in the next step S6, each bending life span corresponding 
25 to maximum stress of each electric wire and/or each wire 
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protecting member, which is obtained in the step S4, is acquired. 
While referring to the predicting functions read out in the 
step S5. From. those bending life spans, a shortest bending 
life span is obtained. The shortest bending life span may be, 
5 for example, the number of bendings for endurance obtained from 
the bending life function in Fig. 4, or may be an endurance 
time calculated based on the number of bendings for endurance 
and the natural frequency corresponding to the former. The 
step S6 corresponds to a bending life predicting step and a 
10 bending life predicting unit in claims. 

in the step S7, a position corresponding to the shortest 
bending life span is specified from its finite element. This 
position specifying is useful in locating a part of the wire 
at Which the wire is possibly disconnected. In a step S8, the 
15 position specified in the step S7 and the bending life span 
predicted in the step S6 are output to the display device 53. 
Because the position on the wire protecting member and/or the 
electric wire corresponding to the shortest bending life span 
is output in addition to the bending life span, more accurate 
20 bending life prediction is ensured. The step S7 corresponds 
to a position specifying step and position specifying means 
in Claims, and the step S8 corresponds an output step and an 
output unit in claims. 

The following procedure may be taken: in a step S9, an 
25 interference part is predicted, and in a step SIO a predicted 
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inference part is output. Specifically, as described in the 
tep S304, the stresses calculated for each natural frequency 
re stored in the result file in a state that those stresses 
are linked with information of the positions (finite elements) 
Of the electric wires and/or wire protecting members. 
Therefore, predictive shapes can be generated for each vibration 
mode from the result file, as indicated by lal in Fig. 7A, la2 
in Fig. 7B, and la3 in Fig. 7C. The position information of 
the Obstacle 2, as described above, is stored in the storing 
device 55. Therefore, if the position information is composed 
with the predictive shape for each vibration mode, one can see 
that there is a possibility that the electric wire will interfere 
with the obstacle 2 (see Fig. 7C) . On the contrary, it is 
possible to predict a part of the electric wire where it will 
interfere with the obstacle 2. The step S9 corresponds to an 
interference part predicting step in claims. 

If in a step Sll, it is judged that a command of outputting 
the result file is present, the contents of the result file 
55c is output in the text format in a step S12. The output 
contents may be displayed by the display device 53 or printed 
on a sheet by the printing device 54. Further, the contents 
to be output may be designated by the input device 52. If the 
contents outputting is not needed, the sequence of processing 
steps may be ended (N in step Sll). 

Thus, the embodiment of the invention succeeds in 



10 



15 



providing a method and apparatus for predicting bending life 
spans of electric wires or the like, which can more accurately 
predict bending life spans of the wires and can satisfactorily 
meet the demand of enhancing the prediction accuracy in an 
environment where vibrations occur, and a program for predicting 
the same in such ways. In particular, a disconnection life 
of the electric wire, which is due to engine vibration of the 
automobile, can be predicted for a short time and without 
performing the endurance test. This is a very advantageous 
effect. Further, also statistically, the bending life span 
is predicted under severe conditions. Accordingly, the 
embodiment contributes to further enhancement of quality and 
wiring path plan. 

It will be readily understood that the invention is not 
limited to the embodiment above mentioned. For example, the 
wire protecting member is not limited to the grommet . Further, 
the invention is not limited to the inside of the automobile, 
but may be applied to a place in a factory where vibrations 
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According to the invention, the plural electric wires, 
the atmosphere temperatures, pre-vibration shapes of the plural 
electricwires, and constraint conditions of the plural electric 
wires are set, and finite element models of plural electric 
wires are formed. Natural frequencies for the pre-vibration 
shapes of the plural electric wires are calculated, and stresses 
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in individual finite elements of the finite element models which 
correspond to the natural frequencies are calculated, and a 
maximum stress is retrieved from the calculated stresses for 
each electric wire . Predicting functions corresponding to the 
plural electric wires and the atmosphere temperatures, which 
are set, are read out. A bending life span corresponding to 
the maximum stress of each electric wire, is acquired while 
referring to the predicting functions read out, and a shortest 
bending life span is obtained from the bending life spans, and 
output. Accordingly, the bending life spans of the plural 
electric wires can be accurately predicted without performing 
the endurance test. As a result, there is provided a bending 
life predicting method which can meet the demands of reducing 
the development period and enhancing the prediction accuracy. 

According to the invention, the plural electric wires 
and wire protecting members, the atmosphere temperatures, 
pre-vibration shapes of the plural electric wires, and 
constraint conditions of the plural electric wires are set, 
and finite element models of plural electric wires are formed. 
Natural frequencies for the pre-vibration shapes of the plural 
electric wires are calculated, and stresses in individual finite 
elements of the finite element models which correspond to the 
natural frequencies are calculated, and a maximum stress is 
retrieved from the calculated stresses for each electric wire. 
Predicting functions corresponding to the plural electric wires 
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and the atmosphere temperatures, which are set, are read out. 
A bending life span corresponding to the maximum stress of each 
electric wire, is acquired while referring to the predicting 
functions read out, and a shortest bending life span is obtained 
from the bending life spans, and output. Accordingly, the 
bending life spans of the plural electric wires inclusive of 
the wire protecting member can be accurately predicted without 
performing the endurance test. As a result, there is provided 
a bending life predicting method which can sufficiently meet 
the demands of reducing the development period. 

According to the invention, the plural electric wires 
are regarded as a wire-like structure of a bundle of the plural 
electric wires, and natural frequencies of the plural electric 
wires are computed. This feature results in reduction of 
calculation amount. Further reduction of the time taken for 
the bending life span prediction is realized. 

According to the invention, when the plural electric wires 
are bundled into a single bundle, and the bundling of the plural 
electric wires is set as one of the constraint conditions in 
the setting step. This ensures more exact prediction of the 
bending life span. Actually, in some cases, the plural electric 
wires are bundled into a single bundle, and in other cases, 
those are not bundled. The bending life span can be predicted, 
while those cases are clearly separated from each other. 

According to the invention, a position specifying step 
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for specifying a position on the wire protecting member or the 
electric wire corresponding to the shortest bending life span, 
is also specified, together with the bending life span, and 
is output. This feature ensures a more accurate prediction 
of the bending life span. 

According to the invention, an interference part induced 
by vibrations on each wire protecting member or each plural 
electric wires is also predicted. Accordingly, with this 
feature, optimum wiring paths can be designed free from wire 
disconnection by its contact. 

According to the invention, a curve representative of 
a lower confidence interval to a population regression function 
statistically calculated using data of the stresses and the 
number of bendings for endurance under" a plurality of typical 
atmosphere temperatures for the wire protecting members and 
the electric wires, is used for the predicting function. 
Accordingly, the bending life span is predicted under severe 
statistical conditions. As a matter of course, a given 
statistical reliability on the predicting function is taken 
into account, and a computing process of the predicting function 
is easy. As a result, thebending life span is severelypredicted 
without addition of complicated processing procedures. This 
contributes to further enhancement of quality and wiring path 
plan. 



